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Abstract. The time-dependent nonlinear kinetic theory for cosmic ray (CR) acceleration in supernova remnants 
(SNRs) is applied studying the properties of the synchrotron emission from SNRs, in particular, the surface 
brightness-diameter (E — D) relation. Detailed numerical calculations are performed for the expected range of the 
relevant physical parameters, namely the ambient density and the supernova explosion energy. The magnetic field 
in SNRs is assumed to be significantly amplified by the efficiently accelerating nuclear CR component. Due to 
the growing number of accelerated CRs the expected SNR luminosity increases during the free expansion phase, 
reaches a peak value at the beginning of the Sedov phase and then decreases again, since in this stage the overall 
CR number remains nearly constant, whereas the effective magnetic field diminishes with time. The theoretically 
predicted brightness-diameter relation in the radio range in the Sedov phase is close to Er (x D~^'^^^. It fits the 
observational data in a very satisfactory way. The observed spread of Er at a given SNR size D is the result of 
the spread of supernova explosion energies and interstellar medium densities. 
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1. Introduction 

Supernova remnants (SNRs) are the main sources of en- 
ergy for the Interstellar Medium (ISM). They also control 
the physical state of the ISM which presumably includes 
the nonthermal component of Interstellar matter, often 
called the Galactic Cosmic Rays (CRs). The synchrotron 
emission of relativistic electrons plays an important role in 
the general study of SNR properties and in CR production 
inside SNRs in particular. All known SNRs are sources of 
radio-synchrotron emission. Several Galactic SNRs were 
recently detected as sources of nonthermal X-ray emission 
which is presumably also of synchrotron origin (e.g. Petre 
et al. 2001). 

The determination of the distances to the Galactic 
SNRs is an important task which is often based on radio 
observations. When there is no direct distance determi- 
nation, estimates can be made by using the radio surface 
brightness-to-diameter relationship (Sr — D). It is how- 
ever not clear whether any functional correlation between 
T,^{t) and D{t) exists for individual objects during their 
evolution in time t, and if so, for which physical reason 
(Green 1984; Case & Bhattacharya 1998). 

CRs are widely accepted to be produced in SNRs 
by the diffusive shock acceleration process at the outer 
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blast wave (see e.g. Drury 1983; Blandford & Eichler 
1987; Berezhko & Krymsky 1988, for reviews). The non- 
linear kinetic theory of diffusive CR acceleration in SNRs 
(Berezhko et al. 1996; Berezhko & Volk 1997) is able to 
make quite definite predictions about the energy spectrum 
and the spatial distribution of CR nuclei and electrons at 
any given evolutionary epoch i, and about the proper- 
ties of the nonthermal radiation produced in SNRs due to 
these accelerated CRs. Application of this theory to indi- 
vidual SNRs (Berezhko et al. 2002, 2003a,b; Volk et al. 
2002) has demonstrated its power in explaining the ob- 
served SNR properties and in predicting new effects like 
magnetic field amplification, leading to the concentration 
of the highest-energy electrons in a very thin shell just 
behind the shock. At the same time there does not yet 
exist a systematic study of the synchrotron emission ex- 
pected during the different evolutionary epochs of SNRs, 
possibly leading in particular to a Er — D relation. An 
interesting question is specifically whether or not the ob- 
served Sr — D relation corresponds to the evolutionary 
track of a single typical SNR, as was for example argued 
by Duric & Seaquist (1986). Such a study is the aim of 
the present paper. 

Previous considerations (e.g. Reynolds & Chevalier 
1981; Duric & Seaquist 1986; Huang et al. 1994) show 
that to explain on average the observed radio synchrotron 
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emission one needs magnetic fields inside SNRs that are 
much higher than those corresponding to typical ISM val- 
ues, -BisM ~ 5 mG. The Raylcigh- Taylor instability of the 
contact discontinuity, separating ejecta and swept-up ISM, 
has been considered as a possible mechanism of magnetic 
field amplification in SNRs (Gull 1973; Fcdorcnko 1983; 
Duric & Seaquist 1986). However, it is unlikely that elec- 
trons which are accelerated at the outer blast wave and 
produce radio emission can penetrate so deeply into the 
interior where the magnetic field is amplified due to this 
instability. The typical energies of such radio electrons is 
so small that their diffusion coefficient is several orders of 
magnitude less than RsVg, the product of shock radius and 
shock speed. This means that these electrons are strongly 
tied to the downstream medium and that their diffusion 
is not important. To play a role in the synchrotron emis- 
sivity of SNRs, magnetic field amplification should occur 
in the same region where also the energetic electrons are 
produced, that is at the shock front. Assuming such a 
scenario Reynolds & Chevalier (1981) concluded that sev- 
eral percent of the shock energy should be converted into 
(amplified) magnetic field energy to explain the observed 
properties of the radio emission during the Sedov phase. 
If the magnetic field was not amplified, the expected time 
dependence of the SNR synchrotron emissivity did not fit 
the observations. 

On the other hand, Luceck & Bell (2000) came to the 
conclusion that a strong magnetic field amplification near 
the shock can indeed be produced nonlinearly by a very 
efficiently accelerated nuclear CR component. In this case 
the CR energy content and streaming anisotropy along 
the field is so high that magnetohydrodynamic waves can 
be strongly excited, amplifying the magnetic field already 
upstream of the shock, where the CR distribution is char- 
acterized by a large spatial gradient. In addition, detailed 
observational evidence has been obtained regarding the 
spatial fine structure of nonthermal X-rays in young SNRs 
(Hwang et al. 2002; Vink & Laming 2003; Long ct al. 2003; 
Bamba ct al. 2003), which can be (explained only if the 
magnetic field at the outer shock and inside SNRs is in- 
deed strongly amplified (Berezhko et al. 2003a; Berezhko 
& Volk 2004). 

In the following we shall present a detailed study of 
the evolutionary properties of the synchrotron emission 
of SNRs. There are a number of distinct differences be- 
tween the present and previous studies: (i) our approach 
is based on the nonlinear kinetic theory self-consistently 
taking into account the backreaction of the accelerating 
CRs on the SN shock structure and remnant dynamics (ii) 
to reproduce the evolution of the synchrotron emission we 
calculate self-consistent electron spectra which take the 
synchrotron losses into account. This is done for every 
SNR evolutionary stage starting from the very early time 
after the SN explosion up to the late Sedov phase when the 
SN shock already becomes an inefficient accelerator (iii) 
we take into account magnetic field amplification near the 
outer SN shock due to CR backreaction on the thermal 
plasma and its magnetic field (iv) to explain the spread of 



the observed properties of the SNR synchrotron emission 
we perform our calculations for different values of the rele- 
vant physical parameters such as SN explosion energy and 
ISM density (v) since recently X-ray synchrotron emission 
was detected for a number of SNRs, we study here the ex- 
pected properties of the SNR synchrotron emission for the 
entire wavelength range, from the radio to the X-ray band. 

2. Model 

A supernova explosion ejects an expanding shell of matter 
with total energy Eg^ and mass Afcj into the surround- 
ing ISM. During an initial period the shell material has a 
broad distribution in velocity v. The fastest part of these 
ejecta is described by a power law 

dM^j/dv oc v"^-^ (1) 

with fc = 7 to 12 (e.g. Jones et al. 1981; Chevalier 1982). 
The interaction of the ejecta with the ISM creates a strong 
shock which accelerates particles. We note that more re- 
cently also exponential profiles for the velocity distribu- 
tion of the ejecta have been discussed (Dwarkadas & 
Chevalier 1998). However we believe that for the global 
synchrotron emission the differences between these differ- 
ent approximations arc not significant. Thus we use here 
the power-law ansatz above. 

Our nonlinear model (Berezhko et al. 1996; Berezhko & 
Volk 1997) is based on a fully time-dependent solution of 
the CR transport equation together with the gas dynamic 
equations in spherical symmetry. Since all relevant equa- 
tions, initial and boundary conditions for this model have 
already been described in detail elsewhere (Berezhko et 
al. 1996; Berezhko & Volk 1997; Berezhko & Ksenofontov 
1999), we do not present them here and only briefly dis- 
cuss the most important aspects below. 

Due to the streaming instability CRs efficiently ex- 
cite large-amplitude magnetic fluctuations upstream of 
the SN shock (e.g. BeU 1978; Blandford & Ostriker 1978; 
McKenzie &: Volk 1982). Since these fluctuations scatter 
CRs extremely strongly, the CR diffusion coefficient is as- 
sumed to be as small as the Bohm limit 

k{p) = K{mc){p/mc), (2) 

where K{mc) = m(?/{2>eB), e and m are the particle 
charge and mass, p denotes the particle momentum, B 
is the magnetic field strength, and c is the speed of light. 

If B is the pre-existing field in the surrounding ISM, 
then the Bohm limit implies that the instability growth 
is restricted by some nonlinear mechanism to the level 
5B '~ B, where 5B is the wave field. The attempt to give 
a nonlinear description of the magnetic field evolution in 
a numerical model ( Luceck & Bell 2000) concluded that a 
considerable field amplification should occur, resulting in 
what we call the effective magnetic field. Broadly speak- 
ing, it was expected that a significant fraction of the shock 
ram pressure PqV^ (po is the ISM density, Vs is the shock 
speed) would be converted into magnetic field energy (Bell 
& Luceck 2001). 
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At the same time, from an analysis of the synchrotron 
spectrum of SN 1006 (Berezhko et aL I2002|l . Cassiopeia A 
(Berezhko et aL 2003b) and Tycho's SNR (Volk et al. 
2002) such a strong magnetic field amplification can only 
be produced as a nonlinear effect by a very efficiently ac- 
celerated nuclear CR component. Its energy density, con- 
sistent with all existing data, is so high that it is able 
to strongly excite magnetohydrodynamic fluctuations, and 
thus to amplify the upstream magnetic field i?isM to an 
effective field Bq > -Bism, and at the same time to per- 
mit efficient CR scattering on all scales, approaching the 
Bohm limit. The same large effective magnetic field is re- 
quired by the comparison of our self-consistent theory with 
the morphology of the observed X-ray synchrotron emis- 
sion, in particular, its spatial fine structure (Berezhko et 
al. 2003a; Berezhko & Volk 2004). 

As the simplest possible conclusion from these inde- 
pendent but convergent considerations we shall assume 
that a constant fraction 5 of the CR pressure Pc at the 
shock is converted into magnetic field energy so that 
B"^ /{%tt) — S ■ Pc- Formally then, the magnetic field in 
the upstream preshock medium becomes time dependent, 
independent of the magnetic field in the ambient interstel- 
lar medium around the SNR: 

Bo = V»7t6-Pc. (3) 

We use the moderate parameter value S — 10^^. In the 
active SNR phase, when P^ ~ Po^si this results in an 
effective downstream magnetic field of energy density 

Bl/{8n) « 10-VoK'- (4) 

Such a large downstream magnetic field value is con- 
sistent with the values derived from the nonthermal X- 
rays measured in individual SNRs (Berezhko et al. 2003a; 
Berezhko & Volk 2004). When in the course of SNR evo- 
lution this field value drops below the magnetic field value 
^ISM in the ambient ISM, we shall use Bq — -Bism = 5 /iC. 
Note that this amplified magnetic field Bq exceeds the 
typical ISM value -Bism ~ 5 /iC during almost the en- 
tire evolution except in the very late SNR evolutionary 
phase. Therefore the results of our calculations below are 
insensitive to the concrete value of Bism- 

To a good approximation the amplified field B^ {r, t) 
in the downstream region is spatially uniform and equal 
to its value B2{t) just behind the shock (sec Appendix). 
This is also roughly true for the free expansion phase. 
For simplicity in what follows we shall make use of this 
property for all evolutionary epochs considered. 

The number of suprathermal protons injected into the 
acceleration process is described by a dimensionless injec- 
tion parameter rj which is a fixed fraction of the number 
of ISM particles entering the shock front. For simplicity 
assume that the injected particles have a velocity equal 
to four times the postshock sound speed. Note that the 
value of the injection velocity Winj is not a significant pa- 
rameter. Since it divides the whole particle distribution 
into two physically different components, thermal particles 



and energetic particles (which we call here cosmic rays), 
the only significant requirement for the value v-mj is that 
it should be high enough, so that the diffusive approach 
is valid for all energetic particles (for a detailed discus- 
sion, see Malkov & Volk 1995). For a given value Ujnj the 
only relevant parameter which determines the number of 
suprathermal particles injected into the acceleration is the 
injection parameter rj. 

Assuming a purely parallel shock, rj is estimated to be 
T] sa 10^2 (g_g^ Scholer et al. 1992; Malkov 1998). As ana- 
lyzed in detail by Volk et al. (2003) the directional struc- 
ture of the fluctuating magnetic field at the shock front 
has the effect of increasingly suppressing the leakage of 
suprathermal particles from the downstream region back 
upstream when the shock becomes instantaneously more 
and more oblique. Applied to the spherical SNR shock in 
a large-scale external field Bism, the local injection rate 
averaged over the fiuctuating magnetic field directions is 
lower than for a purely parallel shock by a factor sa 10"^, 
and even this reduced injection takes place only on some 
fraction /re = 0.15 to 0.25 of the shock surface, depend- 
ing on the size of the SNR. Therefore we adopt here a 
value rj ^ 10^^ for the injection parameter and choose 
a renormalization factor /re = 0.2. The latter is used to 
renormalize (reduce) the number of CRs calculated within 
our spherically symmetric model. These parameter values, 
used for the entire sample of SNRs discussed in this pa- 
per, are close to the values determined individually for the 
three objects SN 1006, Tycho's SNR, and Cas A, men- 
tioned before. In detail the objects in the sample will nat- 
urally be somewhat different. Thus our theoretical analy- 
sis is expected to be appropriate for an entire population, 
even though it will not give a precise result for each indi- 
vidual object separately. 

Since the proton injection rate is not known more pre- 
cisely than by a factor of order unity, we compare the 
results of our calculations performed for rj — 10~^ with 
those corresponding to 77 = 3 x 10^* in order to study the 
sensitivity of the calculations to the value of r/. 

It is assumed that also electrons are injected into the 
acceleration process at the shock front. Formally their in- 
jection momentum is taken to be the same as that of the 
protons. Since the details of the electron injection process 
are poorly known, we chose the electron injection rate such 
that the electron: proton ratio Kcp (which we define as the 
ratio of their distribution functions at all rigidities where 
the protons are already relativistic and the electrons have 
not yet been cooled radiatively) is a constant to be deter- 
mined from the observations. 

The electron dynamics is exactly the same as that for 
protons for electron rigidities corresponding to ultrarel- 
ativistic protons, as long as synchrotron losses are ne- 
glected. Therefore, beyond such rigidities and below the 
loss region the distribution function of accelerated elec- 
trons has the form 

Up) - Kcpfip) (5) 

at any given time. 
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the electron-to-proton ratio, similar to what is observed 

for the Galactic CRs. Only for sufficiently large momenta 
will the electron distribution function /e(p) deviate from 
this relation as a result of synchrotron losses. These are 
taken into account by supplementing the ordinary diffu- 
sive transport equation by a corresponding loss term: 

— = V.V/e-«,V/e + — ^.- + ^-(^-/eJ, (6) 

where the first three terms on the right hand side describe 
diffusion, convection due to the mean gas speed w, and 
adiabatic expansion/compression, respectively. The syn- 
chrotron loss time in the fourth term is given by the ex- 
pression (e.g. Berezinskii et al., 1990) 



(7) 



where rUe is the electron mass and ro is the classical elec- 
tron radius. 

The solution of the transport equations for the ener- 
getic protons and electrons, and of the gas dynamic equa- 
tions at each instant of time yields the CR spectra and 
the spatial distributions of CRs and thermal plasma. This 
makes it possible to calculate the expected nonthermal 
emission produced by CRs in SNRs. 

The choice of K^p allows one to determine the electron 
distribution function and to calculate the associated emis- 
sion. The expected synchrotron SNR luminosity is given 
by the expression (e.g. Berezinskii et al. 1990) 

= 3.8x10-20/,,^ drr'^B^J^ dpp'' U{r,p,t)F (^^^ 

in erg/(s Hz), where 

/>oo 

F{x) =x K5/3{x')dx', (9) 

J X 

Kij_{x) is the modified Bessel function, = 
3eS_Lp^/[47r(mec)^], and B± is the magnetic field 
component perpendicular to the line of sight. It yields 
the flux density S^, at distance d: 



= L,/{47Td^). 



(10) 



Since in the shock region with efficient particle in- 
jection and acceleration the upstream magnetic field di- 
rections can be assumed to be almost completely ran- 
domized due to intense generation of magnetic fluctu- 
ations, we adopt the postshock magnetic field strength 
B2 = crBo, taking into account that at the shock front the 
perpendicular field component undergoes MHD compres- 
sion as a result of the gas compression. Below we shall use 
B± = 0.5B. 



At the start our study, we assume for the SN explosion 
energy a typical value Eg^ = 10^^ erg. We also restrict 
our consideration to the case of a uniform ISM and type 
la SNe which means Mej = IAMq and k = 7 (see Eq. 
(1)). Since the correlation between the SN sites and the 
ambient ISM density structure is not known wc consider 
three different phases of the ISM with hydrogen number 
densities Nu = 3, 0.3 and 0.003 cm~^, which determine 
the ISM mass density as po = l-4?np7VH, where nip is the 
proton mass. The numerical results for the time evolution 
of individual SNRs are shown in Fig.l. 

The gas dynamic part of the problem is characterized 
by the following length, time and velocity scales: 



Ro = (3Mej/4^po)'/', to = Ro/Vo, Vo = ^2E,^/M^;, 

which correspond to the sweep-up radius, the sweep-up 
time and the mean initial speed of the ejecta, respectively. 
The scaling values for the cases considered are i?o = 1-45, 
3.2, 14.9 pc and to = 166, 367, 1709 yr, respectively. 

We note that the following simplified approach is fre- 
quently used for the description of SNR dynamics. It 
treats the ejecta as initially expanding as a whole with 
the single speed Vo- In this case two different SNR evolu- 
tionary phases are distinguished: a free expansion phase 
with roughly constant shock speed Vg Vq (which lasts 
up to t w to) and the subsequent selfsimilar Sedov phase 
for t > to- 

The actual ejecta, created after the SN explosion, have 
a broad distribution in velocity v, cf. Eq. (1). Due to this 
fact the SN shock undergoes substantial deceleration not 
only at t > to but also at t < to- Despite this logical con- 
tradiction we shall still use below the term "free expan- 
sion phase" for the SNR evolution period corresponding 
to t < to- 

As can be seen from Fig.l, the shock speed is never- 
theless constant during a short initial period. This occurs 
because we restrict the velocity distribution of the ejecta 
to a maximum speed Wmax = 4 x 10^ km/s, which is also 
the initial speed of the outer part of the ejecta (piston) 
Vi. This period lasts up to the instant of time when the 
swept-up mass 

= (47r/3)i2,3po (11) 

becomes equal to that part of the mass of the ejecta which 
has the initial speed v = -ymax- During this initial period 
the shock speed is I4 = l.liimax- 

During the subsequent part of the free expansion phase 

{t < to) the shock expansion law is 

Rg OC £;(k-3)/2k^^5-k)/2k^-l/k^(k-3)/k ^^2) 

(Chevalier, 1982) which for k = 7 gives 

Rg cx [El/iM,iPo)]'/h'/'. (13) 

In the following adiabatic (Sedov) phase {t^to) we 
have 



Rs (X (i?.„/po)'/'t'/'. 



(14) 
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t, yr 

Fig. 1. Upstream effective magnetic field Bq, shock radius 
Rs and shock speed Vs (a); total shock (cr) and subshock 
(fJs) compression ratios (b); CR energy Ec as a function 
of time for the case of SN explosion energy Esn = 10^^ erg 
and ISM gas number densities A^h = 3 cm~'^ (solid 
lines), 0.3 cm^"^ (dashed lines) and 0.003 cm^'' (dash- 
dotted lines) in the spherically symmetric case (c). Vi = 
4 X 10^ km/s is the initial piston speed. Thick dashed lines 
correspond to the case Egn — 10^^ erg, A^h — 0.3 cm^'^ 
with higher injection rate 77 = 3 x 10^"*. 



The decelerating shock becomes progressively more 
modified, cf. Fig. lb. The degree of shock modification is 
very well illustrated by the behavior of the shock compres- 
sion ratios. When the CR pressure Pc is small compared 
to the ram pressure paV^ the shock is unmodified and 
has a compression ratio close to the classical value a — 4. 
When Pc becomes comparable to poV^, CR backreaction 
leads to an increase of the total shock compression ratio 
a , whereas the subshock compression ratio as decreases. 



During an initial evolution period t < 10'^ yr the total 
shock compression ratio goes to the value, expected for a 
strongly modified shock, (Berezhko & Ellison 1999) 

a « 1.5M3/«, (15) 

where Ma = Vs/cao is the shock Alfv(e)n Mach num- 
ber, Cao = Bo/y/AirpQ is the Alfv(e)n speed. This limit 
is achieved at the stage t 3 x 10'^ yr. For example, in 
the case of Nh = 3 cm~'^, Ma ~ 30 and ffmax ~ 6, in good 
agreement with the numerical results (Fig. lb). 

Since for t = 10^ - 10^ yr the ampHfied field BqocVs, 
the Alfv(e)n Mach number remains roughly constant. 
Therefore the shock modification is also nearly constant 
during this period of time (see Fig. lb). When Bq even- 
tually becomes equal to -Bism the Alfven Mach number 
Afa oc Vs decreases and, according to Ea. (|15|l . the shock 
compression ratio a oc V^^^ also diminishes in time. 

The non-renormalized CR energy content 

/•OO /'OO 

E,{t)^16TT^ drr^ dpp\i,f{r,p,t), (16) 

JQ Jo 

where ek — ^^Jp^c^ + m'^c'^ — m^c^ is the proton kinetic 

energy, increases during the SNR evolution up to the time 
t ~ 10^ yr, when the shock compression ratio starts to 
decrease rapidly to values cr ^ 4 and even less (Fig Ic) . 

When the SNR age approaches t sa 10*^ yr, the CR 
acceleration efficiency also decreases as a result of the 
fall-off of the Alfven Mach number Ma. Due to its rel- 
atively low speed the shock at this stage continues to 
produce CRs albeit with relatively low cutoff momentum 
Pm(i), whereas previously produced CRs with p > 
leave the shock volume since their diffusive propagation 
is now faster than that of the previously confining shock. 
This process of CR escape from the shock volume can 
be even faster in the case of strong Alfven wave damp- 
ing as a result of their nonlinear interaction (Ptuskin & 
Zirakashvili 2003). Therefore this late stage of the SNR 
evolution, t>3 x 10'* yr, does not play an important role 
in CR production, even though such an old SNR may still 
confine low energy CRs and shine in the radio range. 

4. The accelerated CR distribution 

4.1. Proton distribution 

To illustrate the evolution of the accelerated CR spec- 
trum we consider the case of A'h = 0.3 cm"'^ and Egn = 
10^^ erg. The volume-integrated (or overall) CR distribu- 
tion 

/•OO 

iV(p,0 167rV / drr^f{r,p,t) (17) 

for the CR proton and CR electron components is pre- 
sented in Fig. 2 for different times. The amplitude of the 
spectrum increases with time roughly proportional to the 
CR energy content E^. This process becomes very slow for 
t > 10^ yr as seen from Fig. 2. 
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Fig. 2. The overall CR momentum spectrum for four dif- 
ferent times, for the case of Esn = 10^^ erg and Nn = 
0.3 cm~^. Thick and thin lines correspond to protons and 
electrons respectively. 

The maximum (or cutoff) momentum of the proton 

spectrum at a given evolutionary stage is determined by 
geometrical factors and can be estimated by the expression 
(Berezhko 1996) 

= (18) 
rUpC AK{mpC) 

where 



A = 4[2 + 2b + e- {iy-l)/iy + d]. 



(19) 



The dimcnsionlcss parameters e, b and d determine the 
time variation of the shock size v = din Rg/ dint, of the 
diffusion coefficient e = din k/ dint, of the injection mo- 
mentum b = dlnpinj/dln i?s, and of a measure of the in- 
verse compression ratio d = [21^ — u + Rs{du / dr)]2 / {ui — 
U2). Here u is the plasma speed with respect to the shock 
front and the subscripts 1(2) correspond to the points just 
ahead (behind) the subshock front. 

Let us consider the two SNR evolutionary periods 
when the SN shock is already strongly modified, i.e. the 
end of the free expansion phase and the Sedov phase. For 
a rough estimate one can take Pc = 0.5poV^. In the free 
expansion phase v = 4/7, e = 3/7, b = —3/4, d = —1/2. 
This gives A = 5 and 



Pn 



rripC 



3.6 X 10*^ 



10* km/s 



2/3 



In the Sedov phase we have u = 2/5, e = 3/5, b 
d = 3, which give ^ = 16 and 

4/3 



Pn 



= 1.1 X 10" 



104 km/s 



(20) 
-3/2, 

(21) 



rUpC 

The results presented in Fig. la and Fig. 2 agree very well 
with these estimates. 



Since the most energetic CRs in the power law part 
of their spectrum are produced at the very end of the 

free expansion phase, that is at t ^ tp, the maximum CR 
momentum depends on physical parameters as Pmax oc 
VoRoBoito). This gives 



Pn 



rUpC 



10' 



/ Eg- 



1 



10^1 ergy \1AM(, 



-2/3 



0.3 cm-3 



1/6 



,(22) 



where we have explicitly given the dependence on the 
SN parameters. The maximum CR momentum depends 
strongly on Egn and Mej, but is only weakly dependent 
upon the ISM density pQ. 

We note that a maximum proton momentum Pmax ~ 
lO^TUpC is required to reproduce the spectrum of Galactic 
CRs up to the knee energy (Berezhko & Ksenofontov 
1999). 

The maximum momentum in the overall CR spectrum, 
which is almost constant at the late epoch t > 10*^ yr, is 
about Pmax ~ 10^ nipC and corresponds to the beginning 
of the Sedov phase. All the very high-energy protons pro- 
duced during the preceding free expansion phase t < to 
contribute only a very steep power law tail at larger mo- 
menta p > Pniax in the final overall spectrum. The shape 
of this high energy tail can be derived as follows: 

The kinetic energy of the ejecta contained in the frac- 
tion with speed v within the interval dv is 



dEcj (v) oc 



'dv. 



(23) 



The CR spectrum created by the strongly modified 
shock is quite hard: N{p) oc p~'^, 7 < 2. For the sake of 
a rough estimate wc assume that all the above energy of 
the ejecta goes into CRs with maximum momenta Pm(14), 
where the shock speed is roughly the speed « w of the 
ejecta. This gives 



N{Pm)Pnidpm oc dE^e, 



(24) 



in the ultrarelativistic case. During the period of SNR evo- 
lution when the shock is already modified we have Bq oc Vg 



and therefore Pm oc RgVsB oc Vs 
tain N oc p~'^ with 

7 = (9fc - 33)/6. 



2-!^/(l-!^) 



SO that we ob- 



(25) 



For the case k = 7 this gives 7 = 5 which agrees with 
the numerical results (Fig. 2). This consideration shows the 
following: despite the fact that during the free expansion 
phase the SN shock produces CRs with a power law spec- 
trum up to the momentum Pm(io) which is larger than 
what is attained in the Sedov phase, Pmit > to), these 
particles are seen only in a very steep part of the final 
overall CR spectrum due to their low number (low energy 
content). In addition, significant shock modification does 
not start at the very bcgimiing of SNR evolution but only 
at 100 yr. Due to the fact that efficient CR production 
(followed by strong shock modification and magnetic field 
amplification) takes place only during a relatively short 
period in the free expansion phase, this power law spec- 
trum ends in an exponential cutoff at a momentum which 
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is only a factor of 3 larger than Pmax- Basically these par- 
ticles play no role for the population of Galactic CRs. We 
should add here that in the case of type II and lb SNe, 
which are more numerous in our Galaxy (Tammann et al. 
1994), the value of the parameter k is even larger, /c w 10 
(e.g. Chevalier & Fransson 1984). In these cases the free 
expansion phase presumably contributes even less to the 
final CR spectrum. 

4.2. Electron distribution 

At low momenta we have A^e(p) cx N{p). However, the 
shape of the overall electron spectrum deviates from that 
of the proton spectrum at high momenta p > pi « lO^mpC 
due to the synchrotron losses in the downstream mag- 
netic field which is assumed uniform {B^ = B2 = (jBq). 
According to expression (6) the synchrotron losses become 
important for electron momenta greater than 

The shock continues to produce an electron spectrum 
/c oc p~'^ with g w 4 up to a maximum momentum p^nax 
which is significantly larger than pi. In the momentum 
range from pi to p^jj^x' '^^^ the synchrotron losses the 
overall electron spectrum is therefore close to iVc oc p~^. 

The maximum electron momentum can be estimated 
by equating the synchrotron loss time with the accelera- 
tion time. This results in (e.g. Berezhko et al. 2002) 

Vs 



^ = 6.7 X 10^ 

rUpC 



103 km/s 



(cr-1) /lO^G 



cr(l 



Bo 



(27) 



The main fraction of the electrons with the highest 
energies p ^ lO^mpC is produced at the end of the free 
expansion phase. In this stage, according to Fig. la, Vs ^ 
Vq = 8.5 X 10^ km/s and Bq « 30 ^G, which leads to 
a maximum electron momentum p^^x ~ 4 x lO^WpC, in 
agreement with the numerical results (Fig. 2). 

5. Synchrotron luminosity 

To illustrate the time variation of the synchrotron spec- 
trum during the SNR evolution we present in Fig. 3 the 
calculated synchrotron luminosity L^{v) in an ISM with 
number density A^h = 0.3 cm"^^ for different ages. The 
shape of the spectrum Lv{v) is directly related to the over- 
all electron spectrum Nc{p), shown in Fig. 2. For relatively 
low frequencies < 10^'^ Hz it has a power law form 



(28) 

where the power law index a = (7 — l)/2 is related to the 
power law index 7 which determines the electron spectrum 
Nc ~ AcP~^ near the electron momentum 

p ^ f A^m^c^v V'^ 
m^c \0.87ecB±J ' ^ ' 



which gives the maximum contribution to the synchrotron 
emission at frequency i^. Due to the shock modification the 
electron spectrum has a concave shape characterized by 
an index 7 which slowly decreases with increasing energy. 
Therefore the synchrotron spectrum also has a concave 
shape with a > 0.5 at the lowest frequencies ly ^ 1 GHz, 
decreasing to a = 0.5 at 1 GHz (e.g. Reynolds & 

Ellison 1992; Berezhko et al. 2002). 

The current luminosity L^, is determined by the down- 
stream magnetic field value at the current epoch B^ « 
ctBq , which decreases with time due to the shock decelera- 
tion, and by the total number of electrons produced during 
all the previous stages (x A^. During the initial period, 
t < 10^ yr, the factor A,^ cx A^ grows more rapidly than 
the factor decreases and this leads to an increase of 

in this period. In the later epoch i > 2 x 10^ yr the 
number of accelerated electrons increases only slowly with 
time. Therefore the decrease of the magnetic field strength 
leads to a decrease of L^. 

The synchrotron emission at the highest frequencies 
v > 10^^ Hz is produced by the electrons with p > p\, 
which suffer dramatic synchrotron cooling. Therefore, in 
addition to the above physical factors, the behavior of 
the synchrotron spectrum L^, in this frequency range is 
influenced by the values of pi{t) and Pmax(^)- Taking into 
account expression (|27|) for the electron cutoff momentum 
and relation (|29|l . one can easily derive the value of the 
cutoff frequency 



= 2.1 X 10^^6 



cr-1 



cr2 + 1 ^103 km/s 



Hz, 



(30) 



where b± = B^^j^/B^. It is taken as b± = 0.5 in our cal- 
culation. Since the cutoff frequency is mainly determined 
by the shock speed, its value fmax goes down quickly dur- 
ing the SNR evolution due to the shock deceleration, as 
seen from Fig. 3. It is interesting to note that at frequency 
v ~ 10^^ Hz, corresponding to the typical hard X-ray 
energy = A keV, only relatively young SNRs of age 

103 yr produce very intense nonthermal X-ray emission 
which belongs to the power law part of the synchrotron 
spectrum L^{v). Over the subsequent epochs t'>lQ^ yr 
the shock speed becomes so small that the above X-ray 
frequency range falls into to the exponential tail of the 
spectrum. This leads to a fast decrease of the nonthermal 
X-ray emission during this period of SNR evolution. 

To illustrate the evolution of the synchrotron flux we 
present in Fig. 4 the luminosity Li, dX v ~ \ GHz as a func- 
tion of SNR diameter for the case of SNRs with explosion 
energy £^sn/(10^^ erg) = 0.25, 1 and 3 in an ISM of density 
N^/{\ cm-3) = 0.003, 0.3 and 3. 

The variation of the synchrotron flux is due to varia- 
tion of i?s, Vs and B. 

During a short initial period, which lasts from « 1 yr 
for A^H = 3 cm"3 to t « 10 yr for A^h = 0.003 cm"3 (ggg 
Fig. la), the SN shock speed is constant. Since the SNR 
shock is not significantly modified in the free expansion 
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Fig. 3. Synchrotron spectral luminosity distribution as a 
function of frequency at four different times, for the same 
case as in Fig. 2 (thin lines). Thick lines correspond to the 
injection parameter value r] = 3 x 10~*. 

phase (see Fig. lb), the CR distribution function at the 
shock front has the form 



(31) 



which shows that the CR number density at the shock 
front is proportional to the momentum pinj oc Vg of the 
suprathermal particles injected at the shock front into the 
acceleration process. Therefore we have cx Nb_VsR^, 
Bd oc VIVh and 



L, 



(32) 



taking into account that in the case of an unmodified shock 
7 = 2 and a — 0.5. This explains the rising part of the 
calculated curves Lu{D) which corresponds to small di- 
ameter values D in Fig. 4. 

During the subsequent part of the free expansion phase 
the SN shock decelerates. According to Ea. (|13|l its speed 
as a function of shock radius is 



(33) 



Since the shock, as in the previous epoch, is not yet 
strongly modified, the CR pressure goes like Pc oc NuVs, 
leading to a magnetic field variation oc VTVhK;. The 
expected luminosity dependence on the SNR diameter is 
then 



(34) 



which rather well explains the Li,{D) dependence for D < 
Do = 2i?o. 

In the subsequent Sedov phase a substantial fraction 
of the total energy E^n goes into the CR component whose 
overall number remains nearly constant so that oc _Esn- 
Note that the ISM density iVn influences the amplified 
magnetic field value Bq, which in turn influences the shock 



modification and CR acceleration efficiency (see Fig.l). 
But this influence is not very significant: as one can see 
from Fig.l an increase of TVh by a factor of lO'^ leads to 
an increase of Ec by a factor of only two. Therefore, for 
a rough estimate, E^ and A^ oc E^ can be considered as 
independent of Nn. In this stage the SN shock is signif- 
icantly modified. This is a consequence of the fact that 
the CR pressure is an important fraction of the shock 
ram pressure, with Pc oc poV^. The downstream mag- 
netic field Bd oc u^JTc also varies due to the variation of 
the shock compression ratio a. Since the variation of cr is 
much weaker than that of Po we can neglect it for quali- 
tative analysis, taking B^ oc sfP^ oc \/ A^hKs- Then, taking 
into account the expansion law (|14|l we have 

oc Eli^D-^l^ (35) 

independently of A^h and Afcj, if for simplicity we again 
assume the electron spectrum Ae oc •p~^ with power law 
index 7 = 2, which implies a = 0.5. The calculated elec- 
tron spectra corresponding to the Sedov phase are steeper 
due to the shock modification. This implies 7 > 2 at those 
energies which are responsible for the radio emission. In 
addition, the power law index 7 varies during shock evo- 
lution and therefore the dependencies L^{D) are not of a 
pure power law nature, as is seen from Fig. 4, especially 
in the case of a dense ISM with A^h<;0.1 cm~^ when the 
shock modification is most significant. Nevertheless the 
above analytical relation can be used for rough estimates. 

In the late Sedov phase, when the effective magnetic 
field drops to Bq = Bism and remains constant thereafter, 
the expected SNR luminosity oc Egn is roughly inde- 
pendent of D. 

As pointed out by Case & Bhattacharya (1998), there 
exists no significant correlation between radio luminosity 
and linear diameter. If SNe are more or less uniformly dis- 
tributed over the considered range of ISM densities then it 
is clear from Fig. 4 that no significant correlation between 
Li, and D is to be expected. At least the sample of Case 
& Battacharya appears to be uniformly distributed. 

According to the calculated dependence of on D 
almost all observed SNRs are at the very end of their free 
expansion phase or in the Sedov phase. If this is generally 
so, then the lack of data with small SNR sizes D which 
correspond to the ages t ^ to is a consequence of the 
relatively low expected luminosity of SNRs of age t <^to. 
In addition, it can also be due to a selection effect: SNRs 
with small size and low luminosity are more difficult to 
recognize. 

At high frequencies v ^ 10^^ Hz, which correspond to 
the hard X-ray range, the decrease of the cutoff frequency 
i^iuax becomes appreciable already for t>100 yr. This is 
why the luminosity Li,{D) at 1/ > 10^^ Hz, shown in Fig. 5, 
begins to decrease already in the early Sedov phase. This 
factor becomes dominant for t > 10^ yr and leads to a very 
steep dependence Li,{D). For this reason only relatively 
young SNRs in free expansion or in the early Sedov phase 
are likely to be detected in nonthermal X-rays. As is clear 
from Fig. 5, the relatively high synchrotron luminosity in 
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D, pc 

Fig. 4. The radio spectral luminosity distribution at 
v — 1 GHz as a function of SNR diameter in the 
case of SNR with explosion energy £'sn/(10^^ erg)=0.25 
(thin lines), 1 (normal) and 3 (thick) in an ISM of 
density N-ii/{l cm~^)=3 (solid lines), 0.3 (dashed) and 
0.003 (dash-dotted). Observational data, comprising 37 
Galactic SNRs of known distances, are shown (Case & 
Bhattacharya 1998). 



ent SNR evolutionary stages. When the SN shock is not 
modified (linear CR acceleration regime) , the shape of the 
CR spectrum does not depend on 77 and the amplitude of 
the spectrum is proportional to 77. Such a situation results 
during the initial part of the free expansion phase: as can 
be seen from Fig. lb, at t<10 yr the shock is only slightly 
modified. Since Pc cx; rj, the increase of rj leads to an in- 
crease of the magnetic field value Bq cx (Fig. la), and 
therefore the synchrotron luminosity (Fig. 3) increases by 
a factor of about 7 compared with the case 77 = 10"**. 

When the shock is already strongly modified, i.e. at the 
end of free expansion and in the Sedov phase, the relation 
between 77 and the characteristics of CRs in a SNR become 
more complicated. The shock modification and CR energy 
content, which are closely connected with each other, are 
relatively weakly sensitive to the injection rate as can be 
seen from Fig.l. Similarly, the synchrotron luminosities 
corresponding to the two values of rj considered progres- 
sively approach each other for i>100 yr, so that during 
the Sedov phase t > 10'^ yr the influence of 7/ on Lj^{v) be- 
comes very small (see Fig. 3). Note also that for a higher 
injection rate that leads to larger shock modification the 
synchrotron spectrum in the radio range becomes steeper 
due to the concave shape of the spectrum of accelerated 
CRs. 
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Fig. 5. The X-ray spectral luminosity distribution bX v — 
10^^ GHz (photon energy e,y = 4 keV) as a function of 
SNR diameter in the case of a SNR with explosion energy 
E^^ = 10^1 erg in an ISM of density 7Vh/(1 cm^^) = 3 
(solid line), 0.3 (dashed line) and 0.003 (dash-dotted line). 

X-rays, expected already during the free expansion phase, 
can help to identify very young SNRs that are not yet 
bright radio sources. 

In order to study the sensitivity of our results to 
the injection rate we present in Fig.l and Fig. 3 for the 
case i?sn — 10^^ erg A^h = 0.3 cm^"^ the calculations 
that correspond to the substantially higher injection rate 
77 = 3 X 10""*. The influence of r\ is quite different in differ- 



6. The Yj — D relations in the radio range 

In Fig. 6 we present the surface brightness 

= LJi^'^D'^) (36) 

at frequency v ~ \ GHz as a function of SNR diameter 
D = 2Rs. In the sequel we shall use Er for the brightness 
in the radio band and Sx for the X-ray band, respec- 
tively. Calculations were performed for the above three 
different ISM densities and for the three values of the SN 
explosion energy iJsn/(10^^) erg = 0.25, 1 and 3, respec- 
tively. It was assumed that this energy interval covers the 
majority of the observed SNRs. The experimental values 
corresponding to 37 Galactic shell SNR with known dis- 
tances (Case & Bhattacharya 1998) and to SNRs in the 
Magellanic Clouds, M31 and M 33 (Berkhuijsen 1986) are 
also shown in Fig. 6. 

The relation Tj^{D) is a direct consequence of the 
L^{D) dependence. During the initial part of the free ex- 
pansion phase according to Ea. H32|) the expected depen- 
dence is 

Sr « N^/^D, (37) 

and explains very well the numerical "En^D) behavior for 
small D {D < 0.2 pc for A^h = 3 cm-^ and D < 2 pc for 
Nh — 0.003 cm^^) in Fig. 6. In this region Er does not 
depend on the SN parameters Egn and M^j , but does very 
sensitively depend on the ISM density. 

During the subsequent part of the free expansion phase 
the expected dependence is 

Er oc (£;i/M,j)Vi6Ar2Vi6^-5/i6^ (33) 
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Fig. 6. Surface brightness-diameter diagram at frequency 
v — 1 GHz. Different styles of lines correspond to the same 
cases as in Fig. 4. Experimental data for the Galaxy (Case 
& Bhattacharya 1998) and for the Magellanic Clouds, 
M 31 and M 33 (Berkhuijsen 1986) are shown. 

cf. Eq. H34|l . During this period of SNR evolution the sur- 
face brightness is very sensitive to the values of the explo- 
sion energy Egn and of the ambient gas number density 
Nn, as can also be seen from Fig. 6. 

In the subsequent Sedov phase, according to Eq. H35|l 
we have 

i:nocEU'D-'y\ (39) 

independent of iVfj and Mcj. Since our magnetic field 
strength i3d ex T4 in the Sedov phase changes in a way 
similar to what was suggested by Reynolds & Chevalier 
(1981), the above I]r(Z?) dependence is the same as they 
derived. 

Note again that this dependence Er cx D^^"^^^ is valid 
only for the electron spectrum Nc oc p^^ which is cre- 
ated by the unmodified shock. The actual shock during 
the Sedov phase, as can be seen from Fig. lb, is signifi- 
cantly modified. The modification, characterized by the 
deviation of the shock and subshock compression ratios a 
and CTs from the classical value 4, is larger for denser ISM. 
Therefore the calculated dependence Sr. (I?) presented in 
Fig. 6 is close to Sr (x D~^'^^^ in the case of a diluted 
ISM with number density A^h = 0.003 cm~'^ for which the 



shock is only slightly modified (see Fig. lb) and it becomes 
steeper with the increase of the ISM density. 

In the late Sedov epoch, when due to the shock de- 
celeration the magnetic field amplification becomes irrel- 
evant, and therefore the downstream SNR field becomes 
constant, the Er(£)) dependence goes toward 

Er cx EsnN^/^D-^. (40) 

The fact that the observational data points lie in a 
relatively compact region on the Sr — D diagram can be 
explained from a theoretical point of view if we suggest 
that nearly all of the identified SNRs with known distance 
are in the Sedov phase, or at least at the end of the free 
expansion phase. The lack of a significant number of SNRs 
detected in the early free expansion phase can be due to 
the small SNR size and correspondingly small synchrotron 
flux expected in this stage, as discussed before. 

Even though the expected relation Er cx 
does not depend on the ISM gas density, it is clear from 
Fig. 6 that for small diameters D < 10 pc only SNRs lo- 
cated in a dense ISM with JVh > 1 cm^"^ arc contained in 
the sample, whereas for the largest diameters Z? > 30 pc 
the dominant SNRs are those which are in a tenuous ISM 
with A^H < 0.01 cm~'^. If in addition we assume that SNRs 
in the early Sedov phase are predominantly detected, then 

— 1/3 

one would expect that the relation D cx N-^ is roughly 

17 / 12 

valid, which leads to Sr cx TVjj ' which agrees very 
well with the dependence determined from observations 
(Berkhuijsen 1986). 

Note that type II and lb SNe, whose number is ex- 
pected to exceed the number of type la SNe consider- 
ably, are characterized by essentially different values of the 
mass of ejecta Mcj and of their initial velocity distribution. 
Nevertheless, in the Sedov phase this difference becomes 
irrelevant. Therefore we restrict our consideration here to 
type la SNe. On the other hand, according to expressions 
(|33|l and (|34|l . in the free expansion phase the expected 
radio luminosity of type II and lb SNe with typical ejecta 
masses between Moj = 5 and IOMq, is lower by a factor 
of two to three than that of type la SNe considered here. 

In Figs. 4 and 6 the brightest object is Gas A, which is 
probably a type lb SN with a Wolf-Rayet star progenitor, 
or possibly a type II SN. It is therefore not thought to 
genuinely belong to any of the evolutionary tracks shown. 
However, it remains true that the enormous radio luminos- 
ity and surface brightness of this young and rather small 
object can only be explained by a high-density environ- 
ment, as suggested in these Figs. 

7. The T. — D relation in the X-ray range 

During the last several years a number of SNRs have been 
discovered as sources of nonthermal X-rays (e.g. Petre 
et al. 2001). Studies of the spatial fine structures of the 
nonthermal X-ray emission from SNRs with high angular 
resolution instruments like Chandra (Hwang et al. 2002; 
Vink & Laming 2003; Long et al. 2003; Bamba et al. 
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Fig. 7. The same as in Fig. 6 but for the frequency v = 
10^^ Hz (photon energy ~ 4 keV). Experimental values 
for seven Galactic SNRs (see text) are shown. 



2003) are extremely interesting, since the observed small 
scale structures yield decisive information about the ac- 
tual magnetic field strength in these remnants. The most 
important conclusion is that in all cases the magnetic field 
strength derived from these observations are considerably 
higher than expected for typical ISM fields (Berezhko et 
al. 2003a; Berezhko & Volk 2004). Therefore it is worth- 
while to study the synchrotron luminosity of SNRs in the 
X-ray band. 

In Fig. 7 we show the surface brightness Ex at fre- 
quency 1/ = 10^^ Hz (photon energy ej, = 4 keV) as a func- 
tion of SNR diameter, for the same three values of ISM 
density and SN explosion energy as considered before. The 
experimental brightness calculated from the measured X- 
ray fluxes (and extrapolated to a given energy = 4 keV) 
for Gas A, Tycho's SNR, RCW 86, Kepler and SN 1006 is 
taken from Petre et al. (2001), for G347.3-0.5 from Slane 
et al. (1999), and for G266.2-1.2 from Duncan & Green 
(2000). A distance d = 1 kpc is taken for the last two 
objects. 

In the free expansion phase the T,x{D) dependence is 
very similar to that in the radio range because during this 
epoch the frequency is lower than the cutoff frequency 
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Fig. 8. Spectral surface brightness distribution at fre- 
quency u = 1 GHz versus surface brightness at v = 
10^* Hz (photon energy e,^ Ri 4 keV) for the same cases as 
in Fig.4, except that corresponding to Egn = 10^^ erg. 

t'max (see Fig. 3), and therefore the X-ray luminosity is 
determined by the same factors as in the radio band. 

A dramatically different situation exists for larger di- 
ameters (-D > 4, 8 and 20 pc for A^h = 1, 0.3 and 
0.003 cm~^, respectively), which correspond to the Sedov 
phase. Due to the decrease of the cutoff frequency t'max 
the frequency v = 10^^ Hz becomes larger than i-max al- 
ready at the beginning of this stage. This leads to a very 
rapid decrease of the synchrotron luminosity Li, and of the 
brightness Sx- For example, the brightness Sx becomes 
as low as 10^^ erg/(pc^ s sr) for D = 2 pc in the case 
Esn — 10^^ erg and A^h = 0.3 cm~^ that corresponds to a 
relatively young SNR (see Fig. la) of age t « 2 x 10"* yr. 
It is therefore not surprising that the number of SNRs de- 
tected as sources of nonthermal X-rays is so small, and 
that all of these sources arc young. 

Since the radio and X-ray synchrotron emissivities of 
SNRs are of the same physical nature, it is interesting to 
study the relation to be expected between them. For this 
purpose the radio brightness Y,^ GHz plotted in Fig. 8 
versus the X-ray SNR brightness Sx for the same ranges 
of the explosion energy E^n and ISM number density A'h 
as before. 

It is interesting to note that the allowed region in this 
Sr_Sx diagram is very limited, especially if the ambient 
ISM density is known. This is because a relatively high 
X-ray brightness is expected for only a limited range of 
the SNR sizes, corresponding to the free expansion phase. 
In Fig. 8 all Galactic SNRs except Gas A are inside the ex- 
pected region in the Sr — Sx diagram. A possible explana- 
tion is again that Gas A is not a type la SNR, but certainly 
a type lb or a type II SNR expanding into a nonuniform 
circumstellar medium strongly modified by the progenitor 
star wind. What is even more essential in this case is that 
the mass of the ejecta is substantially higher than in the 
case of a type la SN. This implies a lower shock speed for 
the same explosion energy than in the type la case. As was 
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already pointed out, X-ray synchrotron emission is very 
sensitive to the shock speed. Lowering the shock speed 
will therefore shift the theoretical curve in Fig. 8 mainly 
to the left, toward lower values of Ex- This could lead to 
bounding the experimental Cas A point by the curve that 
corresponds to TVh = 3 cm~^. This mean density appears 
to be an appropriate value for Cas A. 

8. Discussion and conclusions 

The nonlinear kinetic theory for CR production in SNRs 
has been used for the first time in order to study the sta- 
tistical properties of the synchrotron emission of SNRs, in 
particular the Tj — D relation. 

Following the result of theoretical modeling of mag- 
netic field amplification near the shock due to CRs (Luceck 
& Bell 2000), the theoretical analysis of the spectral 
properties of the nonthermal emission of individual SNRs 
(Berezhko et al. 2002, 2003b; Volk et al. 2002) which pre- 
dicted magnetic field strengths inside young SNRs much 
larger than typical ISM values, and the experimental stud- 
ies of the fine structure of nonthermal X-ray emission of 
young SNRs (Hwang et al. 2002; Vink & Laming 2003; 
Long et al. 2003; Bamba et al. 2003), we incorporate into 
our model a time-dependent magnetic field, whose energy 
density is proportional to the CR pressure. This effect in- 
fluences the SNR synchrotron emission drastically, making 
it much more intense, especially in the case of relatively 
young SNRs. 

Since the value of the magnetic field in the most active 
SNR phase for CR production, the early Sedov phase, is 
more than an order of magnitude larger than the ISM 
field, the overall CR proton power law spectrum, produced 
during the whole SNR evolution, extends up to an energy 
Emax ~ 10"'^^ eV, independent of the ISM density. This 
is just the condition required to explain the Galactic CR 
spectrum due to SNRs up to the energy of the knee, which 
is 3 X 10^^ eV (Berezhko & Ksenofontov 1999). 

The maximum energy of CR electrons is about an or- 
der of magnitude lower due to synchrotron losses. In ad- 
dition, as a consequence of these losses, the electron spec- 
trum is much steeper than the proton spectrum at energies 
Ce^lOO GeV. 

Reynolds & Keohane (1999) and Heindrich & Reynolds 
(2001) estimated the maximum electron energy for a 
number of galactic and extragalactic SNRs based on 
the measured radio and X-ray spectra. Their values of 
e^ax ^re within the range of from 10-\/lO jiG/B^ TeV to 
200a/10 fiG/B^ TeV, consistent with our predictions (see 
Fig.2 and Fig. la). 

It was shown that, as a result of the progressively in- 
creasing overall number of accelerated CRs, the radio syn- 
chrotron luminosity increases with time in the free expan- 
sion phase, achieves its peak value at the very beginning 
of the Sedov phase, and then again decreases with time 
due to the decreasing strength of the effective magnetic 
field. The dependence of the radio luminosity on the di- 
ameter, calculated for typical values of the ISM density 



and of the SN explosion energy, covers the region of the 
experimental points on the L,y — D diagram very well. The 
lack of data for small linear sizes £) ^ 10 pc is tentatively 
explained by a low luminosity in the free expansion phase 
that makes it diflacult to recognize SNRs with small size 
and low luminosity. 

In contrast to the radio range, the synchrotron lumi- 
nosity in the X-ray range, which has its peak value at the 
beginning of the Sedov phase as well, goes down much 
more rapidly with increasing D, due to the decrease of 
the electron cutoff momentum. This determines the cor- 
responding value of the cutoff synchrotron frequency. 

The theory predicts a dependence of radio SNR bright- 
nesses on diameter D which is close to Sr = AZ)~^^/* in 
the Sedov phase. All the SNRs corresponding to the same 
explosion energy Esn are expected to align along the same 
fine Sr = AD-^"^/^, independent of ISM density. The am- 
plitude of this dependence scales with the explosion en- 

7 /4 

ergy according to A oc Esn ■ Therefore the variation of 
Esn within a factor of ten, that has to be considered as an 
appropriate assumption for the actual Galactic SNR pop- 
ulation, makes it possible to cover almost the whole region 
in the Sr - D diagram by the lines Er = AD"!'^/*. The 
spread in SNR luminosities at a given SNR diameter D 
is mainly due to the spread of the explosion energy Esn- 
It therefore provides the theoretical basis for the Sr — D 
diagram as an instrument for the distance determination 
to SNRs. 

The expected injection rate r] ~ 10~^ leads to signif- 
icant SN shock modification as a result of the CR back- 
reaction during the Sedov phase. Due to the strong non- 
linear coupling at this stage the final results are relatively 
insensitive to ry. The synchrotron luminosity is sensitive 
to 1] only during the initial part of the free expansion 
phase, when the shock is not modified, scaling roughly 
as (X ff/^. 

It is expected that SNRs in a dense ISM populate the 
part of the line Sr = AD~^'^I^ which corresponds to low 
diameter values D, whereas SNRs in a diluted ISM are 
expected to be on this line for significantly larger D. This 
is an explanation for the empirical correlation Sr oc N^^"^ 
between the radio brightness Er and the ISM number den- 
sity A^H (Berkhuijsen 1986). Since the scale value Rq for 

1 /3 

the SNR size is proportional to A^jj' , the expected cor- 
relation is Sr a N^^^"^. It agrees with the experimental 
relation in a very satisfactory way. 
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Appendix A: 

We consider here the spatial dependence of the effective 
downstream magnetic field in the SNR interior. 
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Since the spatial scales of the amplified magnetic field 
are much smaller than the shock radius Rs, the field can be 
considered random and isotropic. Assuming it to be frozen 
into the thermal gas in the downstream region r < Rg, one 
can write (e.g. Mestel 1965; Chevalier 1974): 

PbP"^/^ = const, (A.l) 

where Pb — B^/{8tt) is the magnetic pressure. Therefore 
at any given radial distance r we have 

PB(r,^)cxPB2(^s)p(r,^)*/^ (A.2) 

where ts is the instant of time when the gas element, at 
the present time t > tg located at radius r, had intersected 
the shock front, r{ts) = Rs{ts). The postshock gas density 
P2 is approximated by a constant value p2 = 4po, where 
po is the preshock density. 

The profile of the downstream mass velocity w{r,t) 
and density p{r, t) of the gas in the Sedov phase can be 
approximated hy w = 3rT4/(4i?s) and p = 4:po{r / Rs{t))^ , 
neglecting here the shock modification by the CR back- 
reaction. The shock radius Rg and the shock velocity Vg 
are given by Rg on t^^^ and Vg{tg) (x tg respectively. 
Therefore the equation of motion of a gas element 

dr/dt^w{r,t) (A.3) 

integrates to 

r/Rs{t) = {tg/tf''^, (A.4) 

which gives tg oc r^^ /t^ . 

Substituting PB2(is) oc Vg(ts) oc tg^^^ and tg{r,t) into 
Eq. HA.2|) . we find that the strong decrease of p"^/^ oc 
into the interior, where fluid elements shocked in the past 
are located, is compensated by the strong increase of Vg 
in the same past, i.e. that B^{r,t) = B2{t). Thus B^{r,t) 
is uniform in this approximation. 
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